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1. INTRODUCTION 

These last years, the renewable sources integration in distribution networks (DN), desired by 
environmental protection demands and the market deregulation, raises the attention for decentralized 
generation problems [1]-[3]. Distributed generation (DG) is routinely called decentralized generation, on-site 
generation, embedded generation, or dispersed generation [4]-[6]. 

Thereafter the insertion of dispersed generation in distribution networks DN, the network voltage 
profile and power losses are affected [7]. The distributed generators contain several kinds of energy sources, 
as combustion gas turbines, induction generators, synchronous generators and renewable energy sources such 
as sun and wind, with modeste quantities. It supplied in the overtaking energy by means of the installation of 
the new distribution network. [8]-[10]. A DG unit should be exploited in a manner that may be exploited in 
parallel with the existing distribution networks [11]. 

The connection of distributed generation is reliant on: [12] 

e the degree of generation previously linked 

e the voltage level where in the DG will be linked 

e the greatness of the DG be linked 

e the category of the DG suggested 

e the robustness of electrical network at the point of connection 

So, DG sizes differed according to the category of users and applications, it can be linked at various 
voltage levels, small sizes are integrated at lowest evels, medium sizes and hight sizes are linked at medium 
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and high voltage levels respectively [12]-[14]. The DG is grouped following the type power supplied to the 

distribution network, in several ranges [15]-[18]: 

-Kind 1 : DG provides only real power: This kind contains microturbines, fuel cells, photovoltaic and batteries. 

-Kind 2 : DG provides real power and reactive power: This kind cotains synchronous machine like 
Combined Heat and Power and gas turbine. 

-Kind 3 : DG provides only reactive power: synchronous compensators are among this kind. 

-Kind 4 : DG provides real power but soak up reactive power: This kind is made up of induction generators 
exploited in wind idustry. 

Previous literature study shows that voltage value of DN should be preserved in a well-defined 
interval, this interval is between the nominal voltage of DN minus 5% and the nominal voltage of DN plus 
5%. It means that the voltage value at all the busses is in this range [19], [20], power losses in transmission 
and distribution lines achieve 4% to 5% of total usage [21]. Thus, the minimization of active and reactive 
power losses in the distribution networks have a great importance for utilities in worldwide [22], [23]. The 
incorporation of dispersed generation unit at unreliable location and size may have adverse effects of 
increased system losses and costs. [24], [25], on other hand optimal placement enhance the DN behavior on 
voltage profile plan, decrease power losses, and increase quality of supplied power [26], [27]. 

Furthermore, several researchers have carried out works on the right site, size and number of DG 
units in DN. These works vary between the analytical and optimization approaches which have perfectly 
determined the ideal site, size and number of DG units. Autor in [28] suggested a method founded on cuckoo 
search algorithm (CSA) to detect the right site and dimension of DG and STATCOM in the radial 
distribution netwoks RDS, the suggested method is tested on several IEEE test systems such as 33 bus and 
136 bus, the obtained results were faced up to other approaches already exist, these results indicate that the 
suggested approach has a precise viewpoint of this crtical issue and causes positive effects on the DN 
behavior. Autor in [29] presented a load flow algorithm based on the backward/forward sweep approach to 
determine load flow issues in RDS, the suggested approach can be used in various DN applications. P. Guru 
et al in [30] used PSO algorithm to investigate the ideal arrangement of DG in IEEE standard 33 bus, the 
obtained outcomes have also been verified by comparing with the previous work discussed in literature 
survey. M. Daneshvaret et al in [31] used exchange market algorithm (EMA) to detect the right DG location 
and size in RDS, EMA applied to solve the DG alloction and size problem on various test systems effectively 
that includes 33 bus and 69 bus IEEE test systems, the numerical results indicated that active and reactive 
power losses are minimized, the voltage profile is improved, and the cost of energy not supplied is 
minimized. K. Adetunji et al in [32] improved the particle swarm optimization (PSO) and the whale 
optimization algorithm (WOA) to have the precise siting and sizing of DG units. Autors in [33] proposed 
gravitational search algorithm (GSA) to find the optimal site of DG and D-STATCOM the algorithm has 
applied on 69 bus IEEE test system and other systems. Another approach adopted by D. Kim and I. Kim [34] 
wherein they proposed a Newton-Raphson algorithm using MATLAB for the optimization of DG units. In 
this paper, we have chosen Newton-Raphson method using Matlab to investigate the impacts of two kinds of 
DG connected at the critical bus on DN behavior. It is analyzed the impact of DGs on voltage profile and 
power losses. 


2. RESEARCH METHOD 
2.1. Load flow analysis 

Power flow study is the significant basic tool to analyse power systems and study their operations [35], 
the aim of load flow analysis (LFA) is to compute the voltage, phase plus real and reactive power of the 
whole busses in the power systems. The buses types in the LFA are swing bus where we know voltage and 
phase, PV bus where we know voltage and real power, PQ bus where we know real and reactive power [36]. 
Newton-Raphson and Fast-Decoupled methods are most known in power flow studies. The nodal power 
system equation is used to drive the basic power flow PF [37]. As shown in (1) is for n-bus system. 


Mil = [Yu] * [v] (1) 


where 

LY; jl is admittance Matrix, 

[4] is currents value at each node Matrix. 
LV; | is buses voltages Matrix. 


Si = Vi * li = Pi + JQ; (2) 
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As results, we have; 


. * * n 
Pi-jQi =V; *=Vi * DG, * V;) (3) 
nb : 
P; = [V; | È> 1%] [Gij;cos( 6; = 6;) + B,jsin( 6; = 6;)] (4) 
= |V;| 2 -W| [Gijsin( 6; — 8) — Bijcos( 6; — 6;)] (5) 
[5], eee , nb 
where 
[V; | : is voltage magnitude per-unit at bus (i). 
TAS : is conjugate current magnitude per-unit at bus (1). 
|P; | : is real power injected per-unit injected at the bus (1). 
1Q:| : is reactive power per-unit injected at the bus (1) 
Gij and Bi; : are conductance and susceptance shunted to the line (ij). 
Ôi : is voltage angle of the bus (1) 
Ôj : is voltage angle of the bus (J) 
Nb : number of buses. 


Power losses in DN vary according to several factors, as the system configuration, the losses level 
across system power lines, transformers, etc. there are two categories of power losses: real power losses and 
reactive power losses. Total power losses in DN are calculated by the (6) and (7) [38]. 


Pas = Sie leer (6) 
Qross = Lier Ui | * xi (7) 
where 
nbr : branches number in system 
Ilil : the current magnitude circulated in branch (1) 


xi andri : reactance and resistance of branch (1). 
Four variables are required to comput power flow parametres, which are P, Q, V and 6 = 6; — ôj. 


2.2. Newton-Raphson method 

It is observable that in the nonlinear (4) and (5), the number of functions and the number of 
variables are unequals, Thus the classical methods cannot solve these equations; which explains the use of 
numerical methods to solve this kind of equations Newton-Raphson method, based on the Jacobian matrix. 
The elements of jacobian matrix are partials derivatives values of P or Q with respect to V or 6, then we write 
the (8) [39]. 
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where: 
AP; and AQ; : are power gap at bus (1). 
Ao; : Voltage angle gap at bus (1). 
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Voltage magnitude and phase angle variations were getting by: 


[45 | | 


[air 


spree] j 
= V [fago] (9) 


Thereby the new voltage magnitude and voltage angle values are: 


yer = VE + Avs 


ôr = bf + Ads 


(10) 
The iterative system stops when the convergence criterion is completed. 
AP | <e 
aQ” <e 
| (11) 


where: 
k is iteration 
E is tolerance 


2.3. Methodology 

An IEEE 57-bus distribution system contains seven generators; forty-two load buses (PQ buses) and 
eighty branches have been considered. The system to study has been operated in Matlab using Newton- 
Raphson method;the IEEE 57-bus distribution system is presented in Figure 1. Lines data and generators 
power limits are provided in Tables 2 and 3 in the appendix [40]. 





Figure 1. Single line diagram IEEE 57 bus test system 
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The methodology includes three stages. The stage number one is to interpret the load flow analysis 
on the power system without DG. The stage number two is to located buses with voltage violation. Avoltage 
violation appears when the bus voltage is less than 95% or more than 105% of per-unit nominal voltage and 
the stage number three is to interpret the load flow analysis with the DG attached at the critical bus. DG 
which provides real power only (kind1) and a DG provides real power and reactive power at the same time 
(Kind2) have been connected separately at the bus with the significant voltage violation (critical bus). 

There are proposed the following scenarios: 

e Case 00: simulation without connecting any DG in the distribution network. 
Case 01: Connecting DG with Ppg = 12.5 MW, Qpe = 0 MVAr and Power factor = 1 
Case 02: Connecting DG with Ppg = 25 MW, Qpe = 0 MVAr and Power factor = 1 
Case 03: Connecting DG with Ppa = 12.5 MW, Qpe = 6 MVAr and Power factor = 0.90 
Case 04: Connecting DG with Ppg = 25 MW, Qpe = 12 MVAr and Power factor = 0.90 


3. RESULTS AND DISCUSSION 

The objectif of this section is to examine the effect of DG connected at the critical bus in distribution 
network. As it was mentioned at methodology sub-section. The load flow analysis was carry out on the 
distribution network without distributed generation, this will be used as a reference case. 

- Case n° 00: The system was operated without linking any DG, the reference case results indicate that 
voltage of busses 31, 33, 46 and 51 listed in Table 1 is not within the 5 percent margin, the rest of 
busses are within the allowed voltage range, the bus number 31 is the bus with significant voltage 
violation, the Figures 2, 3 and 4 illustrate the voltage profile. Real power and reactive power losses are 
equals to 27.86 MW and 121.67 MVAr successively. 


Table 1. Buses with voltage violation in case 00 
Bus Per unit voltage Reference voltage [AV |% 


31 0.936 1.00 6.40 
33 0.948 1.00 5.20 
46 1.060 1.00 6.00 
51 1.052 1.00 5.20 


To evaluate the impact of DG connected at the critical bus on the system power losses and voltage 

profile, the following cases where considered. 

- Case n° 01: a 12.50 MW DG was linked to bus 31, the results indicate that power losses decrease and 
become 25.74 MW and 112.93 MVAr. 

- Case n° 02: a 25 MW DG was connected to bus 31, the results indicate that power losses have continue to 
decrease and become 24.94 MW and 108.20 MVAr. 

- Case n° 03: a 12.50 MW and 6 MVAr DG was connected to bus 31, the results indicate that power losses 
increase little bit and become 25.47 MW and 112.03 MVAr. 

- Case n° 04: a 25 MW and 12 MVAr DG was connected to bus 31, in this case we have got the best results 
24.52 MW and 106.91 MVAr. 


3.1. Voltage profile 
3.1.1. DG kind 01 

The Figure 2 indicate that in the case n° 01 (12.5 MW) the voltage magnitude is improved than in 
the case n°02 (25 MW), increasing DG size is not necessarily a good fact. 


tage 


Vol 











Figure 2. Voltage profile variations with DG kind 01 
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3.1.2. DG kind 02 

The Figure 3 revealed that the voltage magnitude in the case n° 03 (12.5 MW and 6 MVAr) is 
adequate, but in the case n° 04 the voltage magnitude at bus 31 is increased over the international standards 
range (1.128 p u), increasing DG size is not necessarily a good fact in this kind of DG. The Figure 4 also 
show that case n° 03 is perfect than case n° 01. So DG which provides both active and reactive power is more 
useful than DG which provides active power only. 


Figure 3. Voltage profile variations with DG 


Figure 4. Comparison among case n° O1 et case 03 
kind 02 


3.2. Power losses 

Real and reactive power losses without any interconnecting DG are equals to 27.86 MW and 121.67 MVAr. 
It is clear from the Figures 5 (a) and (b) that kind 01 case 01 real power losses were diminished by 7,61 % 
and reactive power losses were diminished by 7,18 %, case 02 real power losses were diminished by 10,48 % 
and reactive power losses were diminished by 11,07 %, kind 02 case 03 real power losses were diminished 
by 8,57 % and reactive power losses were diminished by 7,70 % case 04 real power losses were diminished 
by 11,99 % and reactive power losses were diminished by 12,13 % as DG size is increased as power losses 
are reduced, especially for DG type 02. Therefore, size and kind of a DG play an accurate role in power 
losses. 
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Figure 5. Power losses variations; (a) real power, (b) reactive power 


4. CONCLUSION 

Due the current and future trend of electricity users (suppliers and consumers) is set towards 
increasing the integration of DGs, in particular renewable energy sources, a complete analysis has to be done 
before linking DG units to predict the impacts of DGs kind and size on distribution network behavior, but 
without forgetting the costs. In this paper, the authors used a simple numerical method to compare the impact 
of a DG provides real power and reactive power at the same time and DG provides only real power on 
distribution network voltage profile, real power losses and reactive power losses at the critical bus. Newton- 
Raphson load flow method without and with distributed generation has been used. From the results, it has 
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been shown that integration of the two kinds in distribution networks (critical bus) has allowed enhancing the 
system in terms of power; unfortunately, when we incease the size of DGs a negative effect apreared on 
voltage profile at the critical bus and the neighbor’s buses. Results obtained are confirmed by A.S.O. 
Ogunjuyigbe, et al., [41] and other results available in the literature. Therefore, DG provides active power 
and reactive power at the same time is more useful than DG provides only active power and increasing DG 
size is not necessarily a good fact for distribution networks. 


APPENDIX 
N° F-Bus 
1 1 
2 2 
3 3 
4 4 
5 4 
6 6 
7 6 
8 8 
9 9 
10 9 
11 9 
12 9 
13 13 
14 13 
15 1 
16 1 
17 1 
18 3 
19 4 
20 4 
21 5 
22 7 
23 10 
24 11 
25 12 
26 12 
27 12 
28 14 
29 18 
30 19 
31 21 
32 21 
33 22 
34 23 
35 24 
36 24 
37 24 
38 26 
39 27 
40 28 
41 7 
42 25 
43 30 
44 31 
45 32 
46 34 
47 34 
48 35 
49 36 
50 37 
51 37 
52 36 
53 22 
54 11 
55 41 
56 41 
57 38 
58 15 
59 14 
60 46 
61 47 


Table 2. Lines data 


r 
0.0083 
0.0298 
0.0112 
0.0625 
0.0430 
0.0200 
0.0339 
0.0099 
0.0369 
0.0258 
0.0648 
0.0481 
0.0132 
0.0269 
0.0178 
0.0454 
0.0238 
0.0162 
0.0000 
0.0000 
0.0302 
0.0139 
0.0277 
0.0223 
0.0178 
0.0180 
0.0397 
0.0171 
0.4610 
0.2830 
0.0000 
0.0736 
0.0099 
0.1660 
0.0000 
0.0000 
0.0000 
0.1650 
0.0618 
0.0418 
0.0000 
0.1350 
0.3260 
0.5070 
0.0392 
0.0000 
0.0520 
0.0430 
0.0290 
0.0651 
0.0239 
0.0300 
0.0192 
0.0000 
0.2070 
0.0000 
0.0289 
0.0000 
0.0000 
0.0230 
0.0182 


X 
0.0280 
0.0850 
0.0366 
0.1320 
0.1480 
0.1020 
0.1730 
0.0505 
0.1679 
0.0848 
0.2950 
0.1580 
0.0434 
0.0869 
0.0910 
0.2060 
0.1080 
0.0530 
0.5550 
0.4300 
0.0641 
0.0712 
0.1262 
0.0732 
0.0580 
0.0813 
0.1790 
0.0547 
0.6850 
0.4340 
0.7767 
0.1170 
0.0152 
0.2560 
1.1820 
1.2300 
0.0473 
0.2540 
0.0954 
0.0587 
0.0648 
0.2020 
0.4970 
07550 
0.0360 
0.9530 
0.0780 
0.0537 
0.0366 
0.1009 
0.0379 
0.0466 
0.0295 
0.7490 
0.3520 
0.4120 
0.0585 
0.1042 
0.0735 
0.0680 
0.0233 


b 
0.1290 
0.0818 
0.0380 
0.0258 
0.0348 
0.0276 
0.0470 
0.0548 
0.0440 
0.0218 
0.0772 
0.0406 
0.0110 
0.0230 
0.0988 
0.0546 
0.0286 
0.0544 
0.0000 
0.0000 
0.0124 
0.0194 
0.0328 
0.0188 
0.0604 
0.0216 
0.0476 
0.0148 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0084 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0032 
0.0016 
0.0000 
0.0020 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0020 
0.0000 
0.0000 
0.0032 
0.0000 
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62 48 49 0.0834 0.1290 0.0048 0 
63 49 50 0.0801 0.1280 0.0000 0 
64 50 51 0.1386 0.2200 0.0000 0 
65 10 51 0.0000 0.0712 0.0000 0.930 
66 13 49 0.0000 0.1910 0.0000 0.895 
67 29 52 0.1442 0.1870 0.0000 0 
68 52 53 0.0762 0.0984 0.0000 0 
69 53 54 0.1878 0.2320 0.0000 0 
70 54 55 0.1732 0.2265 0.0000 0 
71 11 43 0.0000 0.1530 0.0000 0.958 
72 44 45 0.0624 0.1242 0.0040 0 
73 40 56 0.0000 1.1950 0.0000 0.958 
74 56 41 0.5530 0.5490 0.0000 0 
75 56 42 0.2125 0.3540 0.0000 0 
76 39 57 0.0000 1.3550 0.0000 0.980 
1 57 56 0.1740 0.2600 0.0000 0 
78 38 49 0.1150 0.1770 0.0030 0 
79 38 48 0.0312 0.0482 0.0000 0 
80 9 55 0.0000 0.1205 0.0000 0.940 

Table 3. Generators power limits 
N° Bus Pmax(MW) Pmin(MW) Qmax(MVAr) Qmin(MVAr) 
1 1 575.88 0 200 -140 
2 2 100 0 50 -17 
3 3 140 0 60 -10 
4 6 100 0 25 -8 
5 8 550 0 200 -140 
6 9 100 0 9 3 
7 12 410 0 155 -150 
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